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Axiomatic principles such as analyticity, unitarity, and crossing symmetry constrain the second 
derivative of the n n scattering amplitudes in some channels to be positive in a region of the Man- 
delstam plane. Since this region lies in the domain of validity of chiral perturbation theory, we 
can use these positivity conditions to bound linear combinations of h and We compare our 
predictions with those derived previously in the literature using similar methods. We compute the 
one-loop TTTT scattering amplitude in the linear sigma model (LSM) using the MS scheme, a result 
hitherto absent in the literature. The LSM values for Zi and h violate the bounds for small values 
of majm-n. We show how this can occur, while still being consistent with the axiomatic principles. 



I. INTRODUCTION 

Quantum chromodynamics (QCD) with two light 
quark flavors has an approximate SU{2)l x SU{2)r chi- 
ral symmetry which is spontaneously broken to its di- 
agonal vector subgroup SU{2)v, leading to an isotriplet 
of pseudo-Goldstone bosons, the pions. Low-energy pion 
dynamics, particularly elastic pion-pion scattering, en- 
codes useful information about the confining dynamics 
of the strong interactions. 

The standard technique to study pion dynamics at very 
low energies with effective field theories was proposed in 
Ref. [H (see also Ref. [H) and systematized as an expan- 
sion in powers of momentum and quark masses in Ref. 3 . 
This effective theory is known as chiral perturbation the- 
ory (xPT). It is formulated in terms of a Lagrangian 
whose only degrees of freedom are pions and which incor- 
porates the symmetries of QCD, including spontaneously 
broken chiral symmetry 0]. 

At lowest order in the chiral expansion, the physical 
observables are determined in terms of two parameters, 
the pion decay constant and the pion mass. If one goes 
beyond the lowest order, a number of low-energy con- 
stants (LECs) li not fixed by symmetries must be in- 
cluded. These can be determined by fitting to experi- 
mental data (for the best determination, see Ref. [j]) or 
estimated by vector- meson dominance but both 

methods have large uncertainties. 

An alternative formulation of tt tt scattering can be ob- 
tained based only on axiomatic principles of quantum 
field theory, such as analyticity, unitarity, and crossing 
symmetry. This allows one to obtain relations between 
observable quantities that must hold, regardless of the 
theory used for the description of the phenomenon un- 
der study. Of course, one of the usual benefits of an ef- 
fective theory approach is that many of these principles 
are automatically satisfied by the scattering amplitudes 
computed using the effective theory. Nevertheless, there 



is still useful information missing in the effective theory 
approach, and one obtains interesting results by study- 
ing the constraints imposed by axiomatic principles on 
the effective Lagrangian. Analyticity and unitarity can 
be exploited to write the well known dispersion relations 
for the scattering amplitudes. These, together with cross- 
ing symmetry, can be converted into positivity conditions 
on scattering amplitudes, which in turn can be combined 
with the xPT predictions to give bounds on the li and 
I2 LECs in the chiral Lagrangian at order p^. 

Two-ffavoryPT was combined with axiomatic princi- 
ples in Ref [7|, which analyzed constraints on s and p 
partial-wave amplitudes in the framework of dispersion 
relations. The analysis was done in xPT at the one-loop 
level. In Ref. Q this study was extended to cover all three 
isospin amplitudes of tt tt scattering at the two-loop level 
in xPT. The best bounds were found for positivity con- 
ditions on full amplitudes (in contrast with partial-wave 
amplitudes), and we follow this approach in the present 
work. However, we find inconsistencies in the domain of 
applicability of the positivity constraints used in Ref. Q 
which will be explained in Sec. [V] Similar bounds were 
first found in Ref. in the context of the Froissart- 
Gribov representation for the scattering lengths. More 
recently, in Ref. [ll|, the very same bounds of Ref. 3| 
were rediscovered using the same procedure as in Ref. 8| 
but using a more restricted domain of validity (in the 
Mandelstam plane) of the positivity constraint. Refer- 
ences d, [Tl| both used one-loop xPT amplitudes. We 
show that the methods of Ref. and Refs. [1, [llj are 
equivalent, and we improve the bounds by properly us- 
ing the domain of validity considered in Ref. [8|, which 
is bigger than that considered in Ref. [ll| . 

In Ref. a different approach was followed for putting 
bounds on some xPT parameters. QCD inequalities on 
Green functions of quark bilinear currents were used to 
obtain relations (inequalities) that involve light quark 
masses, the quark condensate, and some LECs. With 
our method we are insensitive to the quark mass and 
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condensate, since these are lowest order quantities, and 
our analysis starts at C(p^). On the other hand, since 
our study relies on scattering amplitudes, we only make 
use of the chiral Lagrangian when vector, axial-vector, 
and scalar sources are switched off (one always needs the 
scalar source for giving masses to the pions). In fact we 
can only give bounds for the 0{p'^) LECs of operators 
containing only pion fields, li and I2, and so our results 
do not overlap with theirs. 

One of the most popular models used in the litera- 
ture for the study of pion dynamics is the linear sigma 
model (LSM), introduced in the sixties by Gell-Mann 
and Levy [10|. In this model, spontaneous chiral sym- 
metry breaking is driven by a scalar particle a acquiring 
a nonvanishing vacuum expectation value. The LSM La- 
grangian is renormalizable and thus has a reduced (finite) 
number of parameters compared with the most general 
chiral Lagrangian. It shares the same symmetries as xPT 
but has an additional (cr) particle in its spectrum. If the 
(T mass is sufficiently greater than that of the pions, it can 
be formally integrated out of the action, leaving behind 
the xPT Lagrangian, with all the low-energy constants 
having specific values which can be predicted in terms of 
the finite number of parameters of the LSM. 

The values for li and I2 predicted by the LSM do not 
satisfy the dispersion relation bounds for low values of 
the a mass. We will demonstrate that the LSM is per- 
fectly consistent with the dispersion relation bounds and 
that the apparent contradiction results because for low 
values of the a mass, integrating out the a is not valid, 
or equivalently, that higher order terms in the chiral ex- 
pansion cannot be neglected. 

This paper is organized as follows: in Sec. [Ill we de- 
scribe the general features of tt tt scattering such as cross- 
ing and analyticity, and derive the corresponding disper- 
sion relations; in Sec. lIIII we show how dispersive integrals 
imply a positivity condition for the second derivative of 
the scattering amplitude; in Sec. lIVI we convert those pos- 
itivity conditions into bounds for the chiral LECs I12] hi 
Sec. |V| we compare our results with previous analyses, 
and in Sec. IVII we resolve the apparent contradiction be- 
tween the LSM prediction for the chiral LECs and the 
bounds previously found; our conclusions are summa- 
rized in Sec. I VIII In Appendix A we show the relation 
between the methods of Refs. @ and [ll|; in Appendix B 
we calculate the one-loop tt tt scattering amplitude in the 
LSM renormalized in the MS scheme. 



II. DISPERSION RELATIONS FOR tttt 
SCATTERING 

In this section, we find the region of the Mandelstam 
s — t plane in which the tt tt scattering amplitude is 
analytic, and derive the corresponding dispersion rela- 
tions. In a scattering process such as a{pa) + b{pp) 
c{pc)+d{pd) the (nonindcpendcnt) Mandelstam variables 



are defined as 

S = [Pa +Pbf , t ^ {pa- PcY ' U = {pa- Pp)"^ 
4 

s + t + u ^ . (1) 

i=l 

Reversing the order of the two final (or initial) states 
amounts to exchanging t and u. 

We begin by briefly reviewing a few properties of tttt 
scattering. For further details the reader is referred, for 
instance, to Ref. The three pionic states can be 

labeled either by /s = —1, 0, 1 or by Cartesian indices 
a = 1, 2, 3. Both sets of states are linearly related be- 
tween them and to the physical pion states: 

|1,±1) = tK±), |1,0) = |^°), (2) 

where \tt°-) denotes the Cartesian basis, and |l,/3) de- 
notes the isospin basis states. Isospin invariance implies 
that there are only three linearly independent scatter- 
ing amplitudes in the / = 0, 1,2 channels, and crossing 
symmetry relates them to each other, so they can all be 
described by a single function of s and t. In the Cartesian 
basis we can write the Chew-Mandelstam formula 

+ A(u,t,s)rt^% (3) 

where crossing symmetry implies A{x,y,z) ~ 
A(x,z,y) = A(x,y) = A{x,Am? — x — y) where m 
is the pion mass. The function A is related to the isospin 
amplitudes through 

r°(s,t) = ?,A{s,t)+A{t,s)+A{u,s), 

T^{s,t) = A{t,s) - A{u,s) , 

T'^{s,t) = A{t,s)+A{u,s). (4) 

The / = 0, 2 amplitudes are symmetric under the ex- 
change of the final states, whereas the / = 1 is antisym- 
metric: T(0'2)(s^t) =T(o.2)(g^y)^ ri(s,t) = -T1(s,m). 

The isospin amplitudes in the different kinematic chan- 
nels are also linearly related. For our present purposes, 
we only need the relation with the crossed u-channcl. 
This follows directly from Eq. ffl and can be conveniently 
displayed in matrix notation [l3| 

T'is,t) = C'/T''{u,t), 
Cu Cl'' = 61J , 
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(5) 



where, as expected, the crossing-matrix C„ is its own 
inverse. T^{s,t) is the scattering amplitude with isospin 
/ in the s-channcl, and {u,t) is the amplitude with 
isospin /' in the u-channcl. 



3 



Axiomatic principles can be used to show that scatter- 
ing amphtudes are analytic in the full complex s plane 
except for possible isolated points, due to single-particle 
exchange, and branch cuts, due to unitarity. For our 
purposes we only need to know the position sq of the 
first branch point along the real axis of the complex s 
plane. There is then a branch cut along the real s-axis 
for s > So- Any other singularities along the real s-axis 
will be along this cut. The remaining branch cuts will be 
determined by crossing symmetry. 

Let us concentrate on the s-channel keeping t fixed. 
Unitarity ensures that for real s, the scattering ampli- 
tude only develops an imaginary part above the lowest 
mass threshold of possible intermediate states.^ In our 
case the threshold corresponds to two-pion states, i.e. 
So = 4 m^. This means that above the production thresh- 
old (for physical amplitudes) the scattering amplitude is 
complex. Since below threshold, the amplitude is real 
and analytic away from the real axis, it follows from the 
Schwarz reflection principle that T* (s + i e) ~ T {.s — i e) 
and hence T {s + ie) -T {s - ie) = 2ilmT {s + ie) ^ 0. 
This means there must be a branch point at s = 4to^, 
and a discontinuity in the amplitude along the real axis 
for s > 4 771^. We will choose the branch cut to run 
along the real s-axis, because as already explained, the 
other branch points due to higher mass thresholds {e.g. 
four-pion state si = 16 m^) or singularities due to single- 
particle states (e.g. p exchange Sp = will he along 
it. We conclude that our amplitude is nonanalytic for 
s > 4m^, regardless of the value of t. The amphtude 
must also reproduce the singularities in the crossed chan- 
nels, so it is nonanalytic for s,t,u > Arn^. The region 
in the s — t plane where the amplitude is analytic is lim- 
ited to the inside of the triangle defined by the conditions 
s, t < 4m^, s + 1 > 0. Arn? is referred to as the normal 
threshold, associated to the production of two pions. In 
Refs. 0, [m it is assumed that the amplitude is only an- 
alytic between the normal threshold and the abnormal 
threshold, corresponding to s, t, u = 0. The region de- 
limited by the condition < s, i, m < 4m^ is known as 
the Mandclstam triangle (see Fig. [T|) . 

However it has been proved |14l | using very general ar- 
guments that rely on perturbation theory to all orders 
(i.e. that are true for every single Feynman diagram), 
that the amplitude becomes nonanalytic only above the 
normal threshold, and that nothing particular happens 
at s = 0. The region bounded by s,t,u < irri^ is the 
larger triangle shown in Fig. [TJ This is the main differ- 
ence between our method and that of Ref. [ll| . We use 
analyticity in a larger domain, and so obtain more restric- 
tive conditions on the scattering amplitude. Reference Q 
uses the same analyticity domain as we do. However, in 



Above threshold, the physical scattering amplitude is defined as 
the value given by approaching the cut from above, TP'^y''(s, t) = 
T{s + ie,t), with e — > 0. This corresponds to the Feynman ie 
prescription for propagators. 



t 



t-channel 


= 4m^ 


\ 








u-channel 




• 

s-channel 



FIG. 1: Mandelstam plane for tttt scattering. The small 
(blue) triangle in the center is the Mandelstam triangle. The 
big triangle (red and blue area) is the region free from sin- 
gularities. The outer regions (yellow) denote the physical re- 
gions for the three crossed channels. The region bounded by 
the thick black line corresponds to the area A in which the 
positivity conditions are satisfied. 




(a) (6) 

FIG. 2: Contour integrals leading to the fixed-f dispersion 
relations. 

their numeric computations, they include points outside 
this region, which is not justified. 

The derivation of the dispersion relation is quite 
straightforward and is very nicely explained, for instance, 
in Ref. [lH|- For our derivation we consider i as a fixed 
parameter. Wc can then use Cauchy's theorem to write 

T^(s,i) = -l-/dx^^, (6) 

ZTTt X — S 

wherever the amplitude is analytic in a neighborhood (in 
s) of the point (s, t), and where the contour 7 encloses the 
point X — s [see Fig. [^^a)]. Then t < Arri^, and if s > 0, 
we have to use s — > s -I- z e, as already mentioned. From 
the results of Ref. [l3| wc infer that fixed-t dispersion 
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relations hold for t < 4m^, but using solely axiomatic 
principles it can be shown (Ref. [lB|) that they are at 
least valid in the interval — 28to^ < t < 4m^, which 
will be adequate for our purposes. For fixed t, we have 
along the real s-axis a right-hand branch cut for s > 
Aw? and a left-hand branch cut for s < —t. The 7 
contour in Eq. ^ can be deformed into 7', as shown 
in Fig. [^b) in order to express the integral in terms of 
the discontinuity of the amplitude along the real axis. 
In order to do this, the amplitude must fall sufficiently 
rapidly that the contribution from the contour at infinity 
vanishes. If it does not, we can perform n derivatives 
(subtractions) to increase the convergence at infinity. 



d" r 71 ' 

ds" 27ri 



(x-s)"+i ^ ' 



For large enough n that the contour at infinity does not 
contribute, one finds after some straightforward manip- 
ulations that 



ds" ^ ' T I- " 

nil' 
-h(-l)" " 
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{x- s) 



n+1 



{x — u) 



n+1 



ImT^ {x + ie,t). (8) 



The first term is from the discontinuity across the right- 
hand cut. The second term is from the discontinuity 
across the left-hand cut, rewritten using crossing sym- 
metry and Eq. ([5]) to relate the s-channel discontinuity 
in the unphysical region s < to the u-channel disconti- 
nuity in the physical region. 

The best constraint comes from using Eq. ^ with the 
smallest possible value of n. The Froissart bound [l7j 
fixes the minimum number of subtractions needed for 
pion-pion scattering to n = 2. Clearly, if we restrict our- 
selves to s < Am? and s -f- 1 > 0, both denominators in 
Eq. ([8]) are positive, and if n is an even number (for in- 
stance, in our case n = 2) the relative sign is also positive, 
except for the sign of C^^ . 



III. BOUNDS IMPLIED BY THE DISPERSION 
RELATION 

Each isospin amplitude admits a partial-wave expan- 
sion. In the case of spin-zero particles, the amplitude 
depends only on the scattering angle 9, defined as the 
angle between the three-momenta of the first initial and 
final pions, in the center of mass frame. Expanding in 
terms of Legendre polynomials Pi we get: 



T'is,t) = ^(2£-t-l)/,^(s)P,(cos( 



^(2£+l)//(s)P,(l 



2t 



s — Am' 



(9) 



where // (s) denotes the partial-wave amplitudes. The 
optical theorem implies 



Im//(s) = sl3{s)al{s) > 0, 



(10) 



where /3(s) = y 1 — ^-j- is the velocity of the pions in 

the center of mass frame, and aj are the partial- wave 
cross-sections in a given isospin channel. Equation (jlOp 
gives 



2t 



A m' 



(11) 

The partial-wave expansion of the absorptive part con- 
verges in the large Lehmann-Martin ellipse, which, when 
projected onto real s translates into the interval — 4 m^ < 
s < 60 777^. We also need to make sure the absorptive 
part is positive. In Eq. ([5]), the region of integration is 
s > Am'^, and as pointed out in Ref. [3], since Pe{z) > 1 
for z > 1 for all £, if we restrict ourselves to t > 0, 
each partial-wave contribution to the imaginary part is 
positive and so the full imaginary part is itself positive. 
As noted in Ref. (llj . one can find certain linear combi- 
nations J2aiT^ with ai > 0, J2aiCi'^Tj = J2jbjTj 
with bj = J2i ct/ C^l'^ > 0. For these linear combinations, 
the two terms in brackets in Eq. ^ give a positive con- 
tribution. Hence, for these linear combinations, in the 
region A defined as s, t < 4 m'^, t > 0, and s + t > (see 
Fig. [1]) the right-hand side of Eq. fH]) for n = 2 is also 
positive. 

There are three linear combinations which satisfy the 
positivity condition, corresponding to the physical pro- 



7r"'"7r+, and tt+tt'' — > 7r+7r° 



cesses 7r°7r° 

These results are in fact expected and can be deduced 
without any mention of isospin amplitudes. The optical 
theorem ensures that for processes with the same initial 
and final particles a+b a+b, the imaginary part of each 
partial-wave is positive definite. The crossed u-channel 
for those processes has equal initial and final states as 
well, a + b a + b, so for such processes, the imaginary 
part along the right- and left-hand cuts will be always 
positive. The positivity conditions for the three processes 
are 



< 3^T(7r°7r°-.7r°7r°)[(s,0eyl], 
ds^ 

< (7r+7r+ -> 7r+7r+) [{s,t)eA], (12) 



sponding to |r(2)(s,t)+ ir(0)(s,t), ir(2)(.s,i) + 



ir(i)(s,t), and T(2)(s,t), respectively. 
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IV. BOUNDS FOR h AND I2 IN xPT: CHOICE 
OF THE MOST STRINGENT POINT 



It is simple to convert the conditions displayed in 
Eq. into bounds for chiral LECs. The region A cov- 
ers a very low energy domain, and is below the 2 tt thresh- 
old in any of the three crossed channels. In this range of 
energies one expects the chiral expansion to work well, 
so we will approximate the right-hand side of Eq. by 
the xPT result at 0{p^). 

Since the xPT amplitude is derived from a local La- 
grangian, it automatically respects the principles of cross- 
ing symmetry, unitarity, and analyticity. One could 
naively argue that the positivity constraints should also 
be automatically satisfied, but this is not necessarily true. 
As noted in Ref. Q, xPT is an expansion in low mo- 
menta, so the amplitude has polynomial behavior (up 
to logarithms) and grows as or even worse at higher 
orders, violating the Froissart bound. As a result, the 
positivity constraints provide additional information be- 
yond xPT, and give restrictions on the LECs. 

The xPT leading order amplitude is linear in s and 
t and so vanishes on taking the second derivative; the 
next-to- leading order amplitude does not. The 0{p'^) 
amplitude can be found in Ref. 0], and its second 
derivative depends only on two LECs: li and I2 in the 
SU{2)l X SU{2)]i chiral Lagrangian. The amplitude can 
be split into polynomial terms quadratic in momenta 
and masses, and chiral logarithms. The former contain 
the LECs and their second derivatives yield energy inde- 
pendent terms; the latter depend only on momenta and 
masses, are independent of the order LECs, and give 
energy dependent contributions to the second derivative. 
The general structure of the bound can thus be written 



2 

E 

i=l 



a, 



M - fj[is,t)eA] > j = 1,2,3, (13) 



where aji are real coefficients and fj{s,t) are functions 
obtained from chiral logarithms and LEC-independent 
polynomial terms, and j labels each one of the processes 
in Eq. The most stringent restriction is obtained 

for those values of (s, t) that maximize fj{s, t) inside the 
region A: 



Y^a.J, > Mis,t)€A] 



(14) 



It is important to estimate possible corrections to the 
bounds in Eq. (jl4[) coming from 0{p^) terms in the am- 
plitude. The computation of the tt tt scattering amplitude 
at this level of precision was performed in Ref. [131, ^^'^ 
can be split into three pieces: two-loop terms (double chi- 
ral logarithms), that only depend on m and ^V; one- loop 
terms (single chiral logarithms), that depend linearly on 
several 0{p'^) LECs (not only li and ^2); and tree-level 
terms, that depend on 0(/) LECs. In Ref. [1, Eq. H 
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FIG. 3: The h —h region allowed by the positivity conditions 
is shown. The three lines correspond to the three bounds in 
TableP We also show the experimentally fitted values of 
Ref. [jj with their error. 



was calculated with the corresponding 0{p^) amplitude 
for 7r°7r° and tt+tt" at s = 0, t = 4m^. Unfortunately 
the corresponding 0{p^) LECs are badly known (reso- 
nance saturation estimates are usually used), and the 
chiral LECs we want to bound, h and l2, appear again 
in the one-loop terms. In addition the rest of LECs in 
the one-loop terms are symmetry breaking operators, and 
hence appear always multiplied by the pion mass. As a 
result, their numerical values are poorly known. So we 
have only control over the two-loop terms. To get an ed- 
ucated guess for the error from the 0{p^) terms, we will 
multiply the value of the purely two-loop correction by a 
factor of 3. To be more conservative we will adopt as a 
common error for the three bounds, the biggest of these, 
which is 0.4. 

There is one last issue to be discussed before we show 
our results. It is well known that the scalar one-loop two- 
point function is not smooth at threshold (for instance 
its imaginary part is zero below threshold but nonzero 
above). Its first and second derivatives tend to infinity 
when we approach threshold from below. So in order for 
the positivity condition to hold, the coefficients in front of 
these first and second derivatives must always be positive 
below threshold. This is indeed the case in all processes 
under study in our work. 

We find that the maximum of fj(s,t) is always 
achieved for t = Am?, regardless of the process (i.e. for 
i = 1, 2, 3); the value for s does depend on the particular 
process. The maximum of fj is at s = for j ~ 1,2. 
For J = 3, the maximum was found numerically to be 
at s = 1.114 m^. Our results are summarized in Table U 
together with a comparison with the values for the exper- 
imentally fitted LECs h = -0.4 ± 0.6 and h = 4.3 ± 0.1 
from Ref. In Fig. [31 we plot the allowed region in the 
li — 12 space parameter, together with the experimentally 
fitted value. 
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Process 


LECs 


Maximum position 


Bound 


Fit to expt. 


TV TV — > TV TV 
TV+TV'' ^ TV + n° 
TV^TV^ TV'^TV^ 


I1+2I2 
h 

li+3h 


{s = 0,t = 4m^) 
{s = 0,t = 4m^) 
(s = 1.114m^ t = 4m2) 


> = 3.925 ±0.4 
> § = 1.350 ±0.4 
> 5.604 ± 0.4 


8.2 ±0.6 

4.3 ±0.1 
12.5 ±0.7 



TABLE L Bounds obtained by unitarity, crossing, and analyticity and comparison with values extracted from a fit to the 
experimental data given in Ref. j4]. The error on the bound is an estimate of the order p^' terms. 



V. COMPARISON WITH PREVIOUS 
ANALYSES 

As mentioned in Sec. HI there are several studies in 
the literature that combine xPT with axiomatic princi- 
ples. In this section we compare with previous results 
and point out the advantages of the method used here. 

In Ref. 01 only the 7r°7r"_ amplitude was considered, and 
so only bounds on + 2 12 could be obtained. From the 
requirement that the s-wave amplitude has a minimum in 
the interval 1.217 < s/m^ < 1.697 they obtain Ti ± 2 12 > 
3.32 ±0.85. This value is less restrictive than our bound, 
and in addition has a much bigger uncertainty. From 
the once subtracted dispersion relation of the full tt^tt^ 
amplitude they obtain li + 2I2 > 3.3 ± 2.5, which has 
a very large error and is weaker than our bound. Using 
the Froissart-Gribov representation for the c?-wavc partial 
amplitude, they obtained our value for the bound, but 
since a reliable estimate of its error was not found, this 
result was not taken into account in the final results in 
Ref. 0. 

In Ref. 0, the Froissart-Gribov representation for the 
d-wave scattering lengths was used to derive positivity 
conditions. In this way, they obtained the same results 
as us for li + 2I2 and I2, with no errors quoted. In Ap- 
pendix A we demonstrate that this method is equivalent 
to ours for the particular point s = 0, t = 4m^. 

In Ref. Q, the analysis of Ref. Q was repeated, requir- 
ing a minimum of the s-wave amplitude in the same inter- 
val as above, 1.217 < s/to^ < 1.697. Surprisingly Ref. 01 
obtained a much more stringent bound, li +2I2 > 6.16 
(no error quoted) . In view of the discussion in both pa- 
pers, it is our belief that Ref. gives the correct answer. 
The main analysis of Ref. Q uses the same method that 
we do, and in the same domain A. It is argued that the 
most stringent point necessarily lies on the 2 s + i = 4 
line, but we do not see why this should be true. In 
fact, we explicitly find that for the 7r"'"7r"'" amplitude, it 
is not on this line. Furthermore, Ref. 0] only displays 
the bounds at s = (i = Am?), where we get the same 
results for I1 + 2I2 and I2, and at s = — 4 {t = 12 m^), 
where the bounds are much more restrictive. The re- 
sult fi ± I2 > 4.914 quoted in Ref. at s = - 4 
(t = 12m^) is violated by the experimentally fitted val- 
ues of Ref. [J| . Even though Ref. uses the same domain 
A as our analysis, for the numerics, they trespass outside 
this region. The bounds h + 2I2 > 6.923, h > 2.01, 



and h +I2 > 4.914 obtained in Ref. at s = - 4 
(i = 12to^) should not be trusted since the fixed-t dis- 
persion relations are not valid for t > Am?. 

Finally, in Ref. [ll| the same method of Ref. is used, 
but only in the Mandelstam triangle, which is why their 
bound for I1+ZI2 is less restrictive than ours, and does 
not exclude any values for I12 not already excluded by 
the bounds on li + I2 and l2- 



VI. UNITARITY RELATIONS FOR THE 
LINEAR SIGMA MODEL 

In SecHVl we substituted the xPT results into Eq. ^ 
and obtained bounds on some undetermined low-energy 
constants in the effective Lagrangian. One can repeat 
this exercise for theories in which the low-energy effec- 
tive Lagrangian is calculable, to test the validity of the 
bounds. In this section we perform such an analysis for 
the linear sigma model. 

The most straightforward method is to use the predic- 
tions of the LSM for li and I2 for the bounds displayed 
in Table H] As already explained in Sec. HI the LSM is 
invariant under the same symmetries as xPT, and so all 
operators obtained after integrating out the a particle 
must belong to the xPT Lagrangian at some order in the 
chiral expansion. In Ref. (3[ this computation was per- 
formed at the one-loop level and at C(p^), the following 
result was obtained: 



±2 1og^ 

g \ m 
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35 



2 log 



leading to the inequalities 
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6 log 
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191 
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8 log — > 16.94 



(16) 



where g is the (weak) coupling constant of the 0^ term 
in the LSM. It can be written (at leading order) in 
terms of the pion decay constant F-^ through the rela- 
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tion 2g = (m^ — rn^)/ F'^P' These results are obtained in 
weak-coupling perturbation theory to one loop, and have 
corrections of order g from the two-loop graphs. The first 
and third relations of Eq. are always satisfied for a 
weakly coupled theory on which Eq. (fT5|) relies, since the 
24 TT^/g term is larger than the other terms for small val- 
ues g. Note that the coefficient of the 1/ g term must have 
the correct sign for the inequality to be satisfied, which 
it does. The second relation does not involve an inverse 
power of the coupling constant, and is not satisfied for 
large enough values of m/nia- In particular, it is vio- 
lated if rria ^ 4.9 m. One way out of this contradiction is 
that the derivation of the inequality, which relies on the 
Froissart bound, is not valid. But it is not difficult to 
show that the LSM is a local renormalizable theory and 
satisfies the Froissart bound. In the chiral limit to 
and the bound is satisfied. The Goldstone boson is made 
massive by a symmetry breaking term (analogous to an 
external magnetic field). The strength of the symmet- 
ric breaking term must be increased to increase m. The 
symmetry breaking term also contributes to the a mass, 
so another way out is if the region m„/m ^ 4.9 is not 
possible for any values of the parameters in the LSM. By 
explicitly computing the masses in the LSM with a sym- 
metry breaking term, one can show that any jm > \/3 
is allowed, and since -y/S < 4.9 there are allowed values 
for the mass ratio which violate the bound. 

The loophole in the argument is that for low values of 
the a mass, the higher 1/to^ corrections become more 
important. Results in Table |T] rely on the fact that in 
xPT, the scattering amplitude can be safely truncated at 
0{p^), which translates into the statement that the LSM 
amplitude can be truncated at 0{m~^). If is not big 
enough, this approximation receives sizable corrections 
and the chiral expansion breaks down. To violate the 
bound in the second of Eqs. (fT6|) requires to^t ^ 4.9.1 to. 
The chiral expansion is formally an expansion in powers 
of m/rria-, and the bound is violated when /rri^ ^ 0.04, 
a finite distance away from the origin. What is surprising 
is that this number, which is formally of order unity, is 
numerically much smaller than one would have naively 
guessed. 

As a first approach, we include the 1/to* corrections 
to the amplitude and find that then the bounds are vi- 
olated for Too- ^ 5 to, which is not satisfactory, but in- 
dicates that the expansion is slowly converging, 



2 Note that in Ref. Q the relation 2g = m^/F^ is used. Instead, 
we identify i<V with the vacuum expectation value of the cr field 
V, which coincides with the pion decay constant at leading order. 
In the nonlinear parametrization of the LSM, the pion fields are 
collected in the exponential matrix exp(i t-tt/v), which coincides 
with the xPT form after the identification = v is made. In 
addition, m„ depends on the pion mass but at leading order the 
combination — does not. Although in practice the choice 
of Ftt does not affect the results of Eq. I|16|l and the discussion 
in this section, we prefer to use the notation of Ref. [19] . 



and the term contribution moves the result in the 

direction of restoring the validity of the bound, since it 
makes the second derivative of the amplitude less neg- 
ative (see Fig. 2]). To test the LSM bound we will ap- 
ply directly Eq. p^ . rather than the expanded form 
Eq. to the LSM scattering amplitude prediction for 
the TT+TT^ — > tt+tt" process. The second derivative of the 
tree-level amplitude for this process within the LSM van- 
ishes, and so one needs the one-loop result. In Ref. 
this calculation was performed using a mass-dependent 
subtraction scheme. The result is expressed in terms of 
finite two-, three-, and four-point scalar one-loop inte- 
grals, which are then expanded in inverse powers of m^. 
We will use instead the numerical values for the full inte- 
gral expressions. The renormalization procedure followed 
in Ref. [T9| is perfectly acceptable for our computation, 
since the physical amplitudes are scheme independent. 
Most modern computations are done in the MS scheme. 
In Appendix B we give the one-loop LSM amplitudes in 
the mass-independent MS scheme, a result which does 
not appear in the literature. 

Recently, in Ref. [2^, the leading logarithms of the 
scalar-scalar QCD Green function have been calculated 
to two-loop accuracy. For the renormalization the (mod- 
ified) MS scheme is used, as we do in Appendix B. How- 
ever they assume the chiral limit (massless pions) and 
small external momenta, which is equivalent to expand- 
ing in inverse powers of TOo- , approach already followed in 
Ref. [igj . For our analysis we need a nonzero pion mass 
and arbitrary values of the external momenta. Neverthe- 
less, the renormalization program is mass independent 
and so in both calculation should coincide. We agree 
with their results. 

The second derivative d^T(s,4m^)/ds^|s=o for the 
tt+tt" — > tt+tt" process in the LSM is computed for any 
value of the ma/m ratio. The results are shown in Fig.|4l 
which clearly shows that the positivity condition is satis- 
fied at the one-loop level in the LSM for any value of the 
(T mass bigger than the pion mass (even though it would 
suffice to be satisfied for jti^ > \/3m). The apparent 
contradiction of Eq. pB]) was only due to the poor con- 
vergence of the 1 / expansion of the LSM amplitude 
for small m^. The nonlinear sigma model (understood as 
the nonrenormalizable effective field theory obtained by 
integrating the a field out the LSM action) is consistent 
(i.e. obeys the axiomatic bounds) only if we (at least) 
include the 0{p^) contribution. 

This should serve as a warning for the estimate of chiral 
LECs by resonance saturation. In such determinations, 
one starts with a chiral invariant Lagrangian with reso- 
nances as explicit degrees of freedom [5|. The values of 
the chiral LECs are obtained, as in the LSM, by function- 
ally integrating out the hadronic resonances. The ratio 
nip/m ^ 5.5 is of the same order than the value m^/m 
that makes the LSM chiral expansion fail. However, we 
believe that since in the Lagrangians of Q , all LECs are 
already generated at tree level, this anomalous behavior 
is absent. In the LSM, I2 is only generated at one-loop, 
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FIG. 4: Plot of 16 7r^F*d^T(s,4m2)/ds2|^^^ in the linear 

sigma model for the ■k'^-k^ — > -k^tv^ process as a function 
of rric/m. The exact amplitude (blue, continuous line) is 
positive for rria > m. The amplitude up to and including 
terms (red, dashed line), is positive for rrio- > 5m. The 
0{mZ^) amplitude (green, dot-dashed line) remains negative 
for ma < 4.9 m. 

which is why the middle inequahty in Eq. (jl6p does not 
have a 1/5 term and has poor convergence in the 1/™^ 
expansion. 

At this point a natural question arises. Since rriK /m ^ 
3.5 < 5 it could be inconsistent to integrate the kaon 
out of the SU[^) xPT action to obtain the SU{2) chiral 
LECs. In fact using the Li, L2, and La values of Ref. [2^ . 
we obtain h = 5.64 ± 0.84 and I2 == 1.95 ± 0.23 which do 
not agree well with the values quoted in Ref. Q , but are 
in agreement with our bounds. The additional complica- 
tions that arise on imposing the positivity conditions to 
SU{3) xPT are discussed in another paper [23j . 



mer is leading order in the 1/Nc counting and repre- 
sents an expansion in ~ (0.7GeV)~^, the latter is 
subleading in large-iVc and represents an expansion in 
1/A2 ~ (l.lGeV)-2 where ~ 47rF^ and is the 
decay constant of the pion . 

One can use Eq. ([8]) with n = 4 to obtain bounds 
for higher order LECs, using the amplitude up to order 
0{p^). The 0{p'^) LECs in the 0{p'^) amphtude vanish 
on taking the fourth derivative but the one-loop 0{p^) 
chiral logarithmic terms do not. However, one-loop dia- 
grams with one insertion of the C'(p^) LECs contribute 
to terms of order times chiral logarithms, which do not 
vanish on taking the fourth derivative. The 0{p^) LECs 
also contribute to the fourth derivative. Thus one now 
gets inequalities involving the O(p^) LECs and 0{p^) 
LECs plus 0{p'^) chiral logarithms. In addition to hav- 
ing a lot of LECs we no longer compare terms of the same 
order in the chiral expansion. 

The low-energy limit of the linear sigma model La- 
grangian is a theory with spontaneous chiral symmetry 
breaking in which the LECs can be computed in terms of 
the coupling constant g of the LSM. The values of li and 
I2 for this model are in apparent violation of the positiv- 
ity bounds for nia ^ 4.9 to, while the range TOq- > ^/^m 
can be realized in the LSM. We have shown that the ap- 
parent violation is an artifact of the truncation of the 
1/to^ corrections and that the LSM is consistent with 
the positivity conditions for rria > to. 

In a subsequent work [23l | we will apply the same 
method to the SU{?>) case, including the full octet of 
pseudo-Goldstones and generalize the method to take 
into account the explicit breaking of SU{i)v symmetry. 
This gives bounds on Li, L2, and L3. 



VII. CONCLUSIONS 

There are nontrivial constraints which follow from uni- 
tarity, analyticity, and crossing symmetry which must 
be satisfied by any relativistic quantum theory. There 
are some interesting and nontrivial constraints on low- 
energy effective theories which arise by imposing these 
constraints on the effective theory scattering amplitude. 

In this work we have transformed the dispersion re- 
lations for the TT TT scattering amplitude into positivity 
conditions for several processes, valid in a certain region 
of the Mandelstam plane below threshold. This region 
is in fact larger than the Mandelstam triangle, as com- 
monly assumed. These positivity conditions can be con- 
verted into bounds for two LECs of the SU{2) xPT La- 
grangian. Our analysis leads to a stronger bound than 
those obtained previously, since we use positivity in a 
larger region of the Mandelstam plane. The values of the 
LECs extracted from experiment are consistent with the 
bounds derived in this paper. 

One nice feature of the structure of the bounds is 
that it correlates two distinct pieces of the 0{p'^) am- 
plitude: LECs and chiral logarithms. Whereas the for- 
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APPENDIX A: RELATION BETWEEN OUR 
METHOD AND SCATTERING LENGTHS 

In this appendix we wish to demonstrate how the pro- 
cedure followed in Ref. is related to ours. Let us start 
by recalling the definition of the scattering lengths. From 



the partial- wave decomposition of Eq. ([9]), one defines for 
each spin and isospin amplitude the scattering lengths al 



lim 

—►4 



His) 



(Al) 



For even the I ~ 1 scattering length must vanish be- 
cause of Bose symmetry. In Ref. Q these scattering 
lengths can be shown to satisfy the positivity conditions 



2al > 0, 



> 0. 



(A2) 



using the Froissart-Gribov representation. 

It is not difficult to relate the scattering lengths to the 

^-derivative of the total spin-/ scattering amplitude: 

J 4^£! d^r^(4m2,t) 



(2^+1) 



■a 



t=o 



{2e + i) 

where we have used a relation analogous to Eq. ([5]) 



(A3) 
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(A4) 



which follows from crossing symmetry in the i-channel. 

For even £ and / = 1, Eq. (jA3[l implies that the cor- 
responding scattering length is identically zero. To see 
this, recaU that T'^{Am^,t) = -T^{4m'^,-t) by Bose 
symmetry. Now, since the point s = 0, i = 4 lies in 
the region A, for ^ = 2 we know that certain linear com- 
binations of the derivatives appearing in the last equality 
of Eq. (|A3|) must be positive. Inverting Eq. ()A3p we ob- 
tain 



d2 T^{Am^,t) 



di2 



5 

32 



(A5) 



Using the linear combinations that give the amplitudes 
in Eqs. ((T2|), and bearing in mind that 02 = 0, we imme- 
diately reproduce the result shown in Eq. (jA2p plus the 
linearly dependent relation 2 -|- > 0. _ 
We have demonstrated that the method in Ref. [Q] 
corresponds to using positivity at the s = 0, t = A m? 
point in region A of the Mandelstam plane. This is why 
Ref. [6| did not find our third bound, which arises from 
s = 1.114 m^, t = 4m2. 



APPENDIX B: ONE-LOOP nn SCATTERING IN 
THE LINEAR SIGMA MODEL 

In this appendix, we renormalize the linear sigma 
model at one loop for finite pion mass. We will use the 




FIG. 5: One-point a function. Double and dashed lines de- 
note o particles and pions, respectively. The sum of all these 
tadpole graphs must vanish to ensure that perturbation the- 
ory is done around the minimum of the potential, including 
quantum corrections. 



mass-independent MS scheme, instead of the subtraction 
scheme of Ref. [l^ . So our Lagrangian will be split into 
renormalized pieces and counterterms. The basic build- 
ing block containing all fields in the LSM is the SU{2) 
matrix 

Yi = a + i-K-T = v + (j + iTr-T = v + Ti, (Bl) 

where r are the three Pauli matrices, a = a — v is the 
sigma field with zero vacuum expectation value (VEV), 
and V = (a) . The Lagrangian is 



(B2) 



At leading order we get the following relations for masses 
and VEV: 



2 .9 + , 





_ 4P 
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ml 
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ml - 





2ff 



(B3) 



At tree level, one can calculate the VEV directly from 
the Lagrangian by minimizing the potential. At one-loop, 
the most convenient procedure is to impose the condi- 
tion that the one-point g function identically vanishes, 
as shown in Fig. [5l This ensures that we are considering 
quantum excitations around an extremum of the poten- 
tial. It also implies that one-point functions (tadpoles) 
are zero in any graph, so we will not display this topology. 

Let us start calculating the quantum corrections for 
the TT and a propagators, as shown in Fig. [HI The pion 
propagator is diagonal in isospin, and thus proportional 
to (5"^ which we drop. The renormalized one-loop con- 
tributions are 



rj-TJT 



9 



16 7r2 

3g 

16 7r2 



[2{ml-m'')l„^{q^)-2mlA, + 2m''A^] , 
{ml - m2) [I^^iq^) + 3h^{q'-)] , (B4) 
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FIG. 6: One-loop corrections to the pion (top) and a (bottom) propagators. 



where we have defined 



A. = I- log ( — 



[ml A„ - 777,2 ^^j 



1 / mi — m'' mi + m"^ \ , f m'' 

— 1 2 log — 

Z \ s mj. — 777^ / \ 7775. 

V{s) / [S + V{s)f - (7772 - 777^)2 



2s 



7/(s)]2 — (7772 — 7?l2)2 



4 777^ 



4 777^ 



\J[s — (7772 + r772)][s — (777,2 _ 777^)]^ (B5) 

From the rcnormahzation of tlic propagators one can ob- 



tain the running of the g coupling constant 



9 d/i 



3 
2^- 



29- 



(B6) 



which ensures that observables are /i- independent. 

Next we calculate the vertex correction to the a -k -k 
interaction, that is, the irreducible three-point function. 
This correction would affect, among other things, the 
decay of the a into two pions. The diagrams are shown 
m Fig. [71 Since the cr is an isospin singlet, its coupling to 
the pair 7r° tt^ must be proportional to 5°'^ which again 
will not be displayed. The renormalized result then reads 



= -2.97) 



2{ml-m')V„{s) 



-f 6 (7772 - 77i2) 14 (s) + 4 /<,,(r772) 



- 5/7r7r(s) + 3/aa(s) 



where we define three-point one-loop functions as 



(B7) 



1 



2Li2 



2Li2 



47772-27772 +s(/3^(.s)-l) 



87772-27772 +s (3/3^ (s)- 2) 
47772-27772+5(3/3^(5)-!) 



-2Li2 



47772-27772 -I- s (3/3^ (s)-l) 
87772-27772 -1-2 s(/3^(.s)-l) 

-2Li2 



1 



+ Li2 

-I-Li2 



PAml) + (3/?.(s) - 1) iM^l) - 1)_ 

47772-27772+5 (/3^(s) - 1) 
2 7712/3.(7772) + 1^ (/3^(s) - 1) {P,{ml) - 1) 
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2 7^2^.(S) + 1^ (3/3.(s) - 1) (/3.(7772) - 1)J 

/3.(m2) + j^(/3.(s)-l)(/3.(7772)-l)J - 2 ^^2 

' /?^(s)-3 

(/3.(s) - 1) [(3^{ml) + 4^ (3/3,(s) - 1) {I3,{ml) - 1)] 

-27772+5(1-/3,(5)) 1 r -27772+5(1-3/5,(5)) 



(B8) 
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Am'^ -2ml 

s + 2mlPM) - 3s/3,(s) + |^ (3/5,(s) - 1) {PM) + 1), 
(4m2_2m2)s(/?^(.s)-3) 



Li2 
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(B9) 



The last set of diagrams to be considered are the cor- 
rections to the four-pion vertex, that is. the four-point 
irreducible function. Diagrams are shown in Fig [5] and 
only contribute to the tt tt scattering. The structure of 
the amplitude for the process tt" tt'' tt"^ tt'^ is identical 
to that of Eq. ([3]), and our result corresponds to A{s, t, u). 
The renormalized result is 

A^- = -2g + ^\^2{ml ~ m'f [D{s,t) + D{s,u)] 
+ A{ml- m^) [V^{s) + V„{s)] + V„{t) + V^{u) 

+ laa (s) + lh^{s)+2 ^ (t) + /vr (w)] I , 

(BIO) 

where D{s,t) is the scalar four-point one-loop function, 



or scalar box diagram, with all external momenta set to 
to2 and two internal masses equal to m and the other 
two equal to nia-, as can be deduced from Fig. [5] Its ex- 
pression is rather cumbersome and will not be disp layed 
here, but it can be found, for instance, in Ref. [21| . 

All the pieces must be combined together to give the 
one-loop amplitude. First we recall the tree-level ampli- 
tude 



^(s,t)troG-lovol = -25 J, (Bll) 

s — 

which reduces to the well-known 0{p^) xPT result when 
the TTio- —I- oo limit is taken. The renormalized one-loop 
amplitude is then 



Ms, t)i_ioop = i 2 {ml - m^y [D{s, t) + D{s, u)] + 4 {ml - m^) [V^{s) + V,{s)] + V,{t) + V,{u) 



2 



8 7r2 



3 (mi — m^y 

-Iaa{s) + 7/^^(s) + 2 [I^-,{t)+I^^{u)] + -^-^ [h r.{s) + i Iaa{s) - 3 A^] 

[ml - s) 

2 {ml ~ m^) Va{s) + 6 (m^ - m^) V^{s) +Ah„{m^) + 5^^(s) + 3Iaa{s)] \ , 



(B12) 

and the total amplitude to one loop is given by adding the two. It is /i-independent once we take into account the 
running coupling constant of Eq. (|B6[) . 
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